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Abstract 

Using the extension of the standard Hawking-Hartle prescription for defining a wave function 
for the universe, we show that it is possible, given a suitable form for the scalar field potential, to 
have the universe begin at its largest size and thereafter contract, with the growth of perturbations 
proceeding from small, at the largest size, to largest when the universe is small. This solution 
would dominate the wave function by an exponentially large amount if one chooses the Hartle 
Hawking prescription for the wave-function, but is exponentially sub-dominant for the Linde- 
Vilenkin prescription. 
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I. LARGE INITIAL UNIVERSE 

In a surprising result, Hartle and Hawking twenty years ago showed [|TJ how, at least in the 
semi-classical theory of quantum gravity, one could give an intuitively appealing definition for 
a preferred wave function for the universe. If one believes that nature selects this definition, 
then they argued that this wave function explains various features which we observe the 
universe to have. One of these features is that the universe is simple when it is small || and 
becomes more and more complex as it grows in size. While Hawking originally believed that 
the universe must always be simple when small || , Page and Laflamme [|J pointed out that 
this wave-function would be expected to also contain probablilities for the universe to be 
complex when small. This was interpreted as indicating a classical history of the universe 
which began simple when small, expanded, with the gravitational instability causing growth 
in the quantum fluctuations present when small. The universe would eventually re-collapse 
with the fluctuations continuing to grow. If one wrote the HH wave-function in terms of 
a basis which contained such semi-classical components, one would interpret them to say 
that the thermodynamic arrow of time would be driven by these fluctuations and that one 
should interpret those solutions such that the low fluctuation end should be interpreted as 
the beginning and the high as the end of the evolution. 

However, there is an alternative interpretation of the HH wave-function ( see for example 
Jheeta and Unruh[||). The prescription bears a strong resemblance to the prescription of a 
semi-classical WKB wave-function near a classical turning point of a potential. The classical 
solution for the HH models is that of a De Sitter space with the the cosmological constant 
dominating the evolution- ie the classical solution looks like De Sitter space in spherical 
coordinates near the point where the spherical space reaches it minimum size. 

Interpreted in this way there is another place where the universe could and does bounce; 
namely if a closed universe has matter such that the energy density of the matter drops as 
the universe expands. In this case, the curvature terms in the FRW equations eventually 
dominates and halts the expansion and the universe re-collapses. The question which we 
address is whether the HH prescription could be used to place boundary conditions at this 
bounce point rather than at the point where the universe achieves it minimum size? 

We will examine a simple model, with a homogeneous scalar field driving the dynamics 
of the universe. The space-time will be assumed to be a closed spherical homogeneous and 
isotropic FRW cosmology with metric 
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The mini-superspace action for this space-time becomes 





which has as semi-classical equations of motion 
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where ' — (The third is the constraint equation, and is an integral of motion of the other 
two). 

The semi-classical HH prescription is that the path integral will be dominated by a 
solution of the equations which has only a final boundary given by the classical state on 
which one wishes to evaluate the wave function ,(I.e., the single boundary is defined by the 
values of a and (or 0) which one wants on the final boundary). We must find solutions 
to these equations which have no initial boundary- which have a(t) going to zero and 0(t) 
regular as a(t) approaches zero. This is impossible if iV ( and a) are real. We can however 
find such solutions if N is chosen to be complex, and in particular if we choose N to become 
purely imaginary as a approaches zero. (For a generalisation to arbitrary complex N see 
Jheeta and Unruh ||). One choice, if it works, is to take N and a to be purely real over 
some range of times (the Lorentzian regime), and then where a(t) has zero derivative, to 
take iV to be purely imaginary. This allows one to smoothly join a(t) across the surface 
where N discontinuously. 

However, in general, both a and cannot be taken to be zero at the same time t and still 
find a solution with only a single boundary. Since the equations of motion for a and are 
smooth second order equations, the derivatives j^r must also be smooth. While we can keep 
a(t) real across this boundary, since going from real zero to imaginary zero is a continuous 
transformation, we cannot do the same for 0. Thus, must thus become complex as one 
crosses this boundary. 

We are going to try to model the universe such that at time t = we have a = 0, and for 
t > 0, the universe contracts from the value of a at this time. Ie, we take t = to be the 
time at which the universe achieves its maximum value in the Lorentzian regime. 

Examining the constraint equation for a, let us examine this boundary where a(t) = 0. 
Since a is, by assumption less, than zero (the universe is contracting) for larger values of t, 
the potential term must be sub-dominant. The curvature term 6a cancels with the kinetic 
energy. If V were to dominate, then a 2 V decreases as a decreases, and a 2 would become 
imaginary. Thus, in order that we have such a maximum in the radius of the universe, the 
term a 3 2 must be dominating the energy, and 2 must be decreasing at least as fast as 1/a 2 
But for very small V this is precisely how is doing. (For small obeys 

T<*£r - (6) 

dt Ndt y ' 

or 

— - 1 (7) 

Ndt a 3 1 ; 

Ie, in order that the universe re-collapse, the potential term cannot dominate the evolution 
of the universe- it must be dominated by the kinetic energy of the matter. 

Now, -k<j) must be continuous across that boundary where N goes from real to imaginary. 
Since we are taking N to be purely imaginary for t < 0, must also be imaginary for times 
just before zero. Ie, must pick up an imaginary part. 

Let us now examine the evolution of this model in this realm where iV is purely imaginary. 
Let us assume that the potential V(0) is such that the potential remains real for = o + i$, 
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where 0o is the value of at the point where d = 0. Since N is arbitrary, we will take N — 1 
for t > and N = i for t < 0. The equation of motion for $ will be 




Tr a ~dt = rf$ (8) 
and for a, the constraint equation gives 

- 1 + ^f 2 + ^*» < 9 > 

The second order Einstein equation for a shows that through t — 0, a reverses sign. Thus, 
while the universe contracts toward the future in the Lorentzian regime, it must expand into 
the past at this turning point. 

Thus ^ is less than zero for t < near 0, and the universe expands as go to negative 
t. This raises the question of how one could get a to go to zero in this N-imaginary region. 
The kinetic energy which dominates is driving the universe larger. The answer is that as 
a increases, $ keeps decreasing. If the potential V{<&) remains non-zero, eventually it will 
dominate once again. Once it dominates, the potential term will keep increasing until it 
becomes of the same order as the curvature term (the 1 in the equation for d) and a has 
another turning point. The universe will then collapse as a function of t, the kinetic term 
will again grow, and the danger is that one will hit another turning point. If however 
the potential is arranged so that during this collapse, the potential increases, and thus $ 
decreases, one can arrange the system so that a can collapse all the way to zero smoothly. 

Let us examine a specific model. We take the potential to be 

V(<i>) = V + a<p\ (10) 

We are going to assume that at t=0, 0(0) = O = 0, a(0) is large, and d(0) = 0. We will 
assume that Vo is smaller than l/a(0) 2 and thus <fi 2 ~ -^yz- 

If the solution is such that a approaches zero smoothly at t — to (so that $ and $ are 
both finite there) then near t — t a we have 

$(t)^$(t )(l--(t-to) 2 + ...) (11) 
a{t) = (t- t )(l - (V - !<&(f ) 4 )(* - t ) 2 + ...) (12) 
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Figure 1 

The evolution of (p and a in time in the N = i regime. We have Vq = 10~8 ; a the quartic 
coefficient is 2 1CT 11 . The value of a at t = is carefully tuned. 

In figure 1 we plot a solution with V = 1CT 8 and a = 2 1CT 11 . The lower axis is Note 
that the field $ approaches a constant for most of the evolution. 

If the potential in the imaginary direction is sufficiently flat for a large enough value of 
$, the universe can oscillate between the the kinetic $ dominated bounce at smaller values 
of a and the potential dominated one at large values of a for many bounces . In Figure 2 
we have a solution with V = 1CT8 and a = 10~ 13 , where we have three bounces before the 
radius eventually goes to zero. 

In the N real regime, the potential a will decrease, and the kinetic dominated term will 
grow roughly as = 1/a 3 . If the kinetic energy dominates, the solution for a goes as t 1 /3, 
and grows logarithmically. If the potential in the real direction grows, it may eventually 
dominate the matter again, and the universe would suffer a De Sitter like bounce at small 
values of a. 

Figure 2 is the solution for N and real into the past from the large radius bounce 
which matches to the solution for Figure 1. In this particular case, because the po- 
tential eventually grows so rapidly with 0, the smaller bounce radius is not that much 
smaller than the largest size. However this depends entirely on the behaviour of the po- 
tential for large real 0. Thus, if our potential were say V cos(0(£)/0 a ), with (f> a very 
large, our potential along the real axis would remain bounded by O , and the dynam- 
ics of a would be kinetic dominated for all times to the past of the large bounce. In 
the imaginary direction, V(<f>) would have the required confining shape and would, once 
had grown to of order a , allow for the radius to go smoothly to zero. In the real 
and iV regime, the solution would be kinetic dominated all the way to a singularity. 
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Figure 2 

Here alpha is taken as 10~ 13 , with the rest of the parameters the same as in figure 1. The 
initial value of is carefully tuned to give the multiple oscillations for a 

In the spirit of HH, we would interpret this universe as beginning in the Euclidean regime 
with a = and <fi = i$ and regular. The universe would "expand" and bounce between 
the potential and kinetic dominated bounce points for as long as desired, until it suddenly 
makes the transition to real N and real at the lower kinetic dominated bounce point. The 
universe would then collapse into the future, until finally, either the potential dominated and 
the universe suffers a De Sitter, inflationary type bounce, or until the the universe finally 
collapses to a singularity. 



II. SIGN OF EUCLIDEAN ACTION 

Since these solutions are possible, how would they enter into the path integral for the 
"wave function of the universe"? There has of course been much heat expended on the 
various schools as to the sign one should take for the action in the Euclidean regime. These 
amount to either a choice of sign for N in the imaginary regime, or a choice in the sign 
of a (since a enters only quadratically in all expressions of physical relevance). Thus the 
contributions of these imaginary N parts of the action will contribute either e s or e~ s where 
S is the integrated action for the imaginary N regime (or is the imaginary part of the action 
in the case of a generic complex metric and field dominating the action). For our purposes, 
this will be the value of 

S I = J (-6a + ad 2 + a 3 Q$ 2 + V(I$)^ alt (13) 

During the bouncing, all the terms are of the same order as the first term, and we have 

S T oc J adt (14) 

By hypothesis both a and t are large, and furthermore, t in general will be of the same order 
as na, where n is the number of bounces in the imaginary iV regime, making Si of order of 
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na 2 . a here will be of order its maximum where a 2 oc y^. Ie, the action will be of order For 
small Vq, and for the positive choice of the S in the action, this would completely dominate 
the path integral. The theory would predict with overwhelming probability that the universe 
began large and then collapsed down to a small size, rather than the more traditional way 
around. Furthermore, the fluctuations would grow during the collapse phase. The universe 
would grow more and more chaotic as its size decreased, with the simplest being when the 
universe was at its largest diameter. This sign is the sign usually attributed to Hartle and 
Hawking [|TJ. 

The alternative sign, usually attributed to Linde and Vilenkin || , J7j] (although derivable 
as simply one of two possibilities in defining the Hartle Hawking wave function) would on the 
other hand make the contributions of this particular solution to the semi-classical equations, 
completely sub-dominant. 

One of the interesting features of this analysis is the realization that the in order to 
maintain continuity of the scalar field equation through the transition from the real to 
imaginary N, one must also allow (f> to become complex. In particular, if the field and scale 
factor are real in the real N regime, then the derivative of the field in general must be 
pure imaginary in the imaginary N regime. Thus, we note that the sign of the kinetic term 
0,3 2W0 2 does not change on the transition from real to imaginary <fi. 



III. FLUCTUATIONS 

Let us examine the action for the inhomogeneous minimally coupled scalar field ip in the 
geometry. The action for this field is 

,3 ( 1 J,2 1 (^l\2 \ A +AC\l 



S = J Na - ^a(WJ dtdn ( 15 ) 

where the dfl 3 is the spatial integral, and (d<p) 2 represents the spatial derivative part of the 
action. Again at the surface where N makes the transition from real to imaginary, and a is 
zero, ip must again be complex to keep ip/N continuous. Ie, again we can write ifi = ipQ + i^f. 
The action then becomes, taking N = i, 

S e = i j a\U> 2 + ^(d^) 2 )dtdn 3 (16) 

Taking the semi-classical approximation to the path integral with the VL sign we have e lSe 
as the contribution to the wave-function. All the terms in the action are positive, indicating 
that we want to minimize the action to find the dominant terms. Using the apherical 
eigenstates for the spatial modes of the field we have the wave function factor as 



e 



-f^^n+^nw (17) 



near a = 0, a will again be linear in t, and the solution for the field which is regular near 
t = t a is 

^^ (l + ^ 2 + ... (18) 

o 

The action exponent can be written as 

/ a3{ ¥ 1 + ^ l)dt = « 3 **(*)**(f)ILo ( 19 ) 
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k 2 



The lower limit is zero. As long as > ^ will "slow roll" and 

q k (20) 



3ad 

while in the other limit, 

* k « (-)** (21) 
a 

Thus, at the transition point, the exponent will be proportional to — ty\ (ie quadratic) and 
will be concentrated near zero with a width proportional to either for short wavelengths, 

or y -^- for long wavelengths. 

Thus, for the large initial conditions, the universe will start in its vacuum state at the 
largest turning point of the theory. 
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